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Abstract 



This paper is devoted to a statistical analysis of the fluctuations of veloc- 
ity and acceleration produced by a random distribution of point vortices in 
two-dimensional turbulence. We show that the velocity probability density 
function (p.d.f.) behaves in a manner which is intermediate between Gaussian 
and Levy laws while the distribution of accelerations is governed by a Cauchy 
law. Our study accounts properly for a spectrum of circulations among the 
vortices. In the case of real vortices (with a finite core) we show analytically 
that the distribution of accelerations makes a smooth transition from Cauchy 
(for small fluctuations) to Gaussian (for large fluctuations) passing probably 
through an exponential tail. We introduce a function T(V) which gives the 
typical duration of a velocity fluctuation V; we show that T(V) behaves like 
V and V^ 1 for weak and large velocities respectively. These results have a 
simple physical interpretation in the nearest neighbor approximation and in 
Smoluchowski (1916) theory concerning the persistence of fluctuations. We 
discuss the analogies with respect to the fluctuations of the gravitational field 
in stellar systems. As an application of these results, we determine an approx- 
imate expression for the diffusion coefficient of point vortices. When applied 
to the context of freely decaying two-dimensional turbulence, the diffusion 
becomes anomalous and we establish a relationship v = 1 + | between the 
exponent of anomalous diffusion v and the exponent £ which characterizes 
the decay of the vortex density. This relation is in good agreement with 
laboratory experiments and numerical simulations. 



PACS numbers: 47.10.+g, 47.27.-i 
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I. INTRODUCTION 



A basic problem in fluid turbulence is the characterization of the entire stochastic varia- 
tion of the velocity field V(r, t) produced by the disordered motion of the flow. The velocity 
fluctuations can be described by different quantities such as their probability density func- 
tion (p.d.f.), their typical duration and their spatial or temporal correlations. We consider a 
simple model of two-dimensional turbulence for which it is possible to calculate these quanti- 
ties exactly. In our model, the velocity is produced by a collection of point vortices randomly 
distributed in the domain with uniform probability. Point vortices behave like particles in 
interaction and share some common features with electric charges or stars []. In particular, 
the problem at hand is directly connected with the problem of the fluctuations of the electric 
field in a gas composed of simple ions or the fluctuations of the gravitational field produced 
by a random distribution of stars. These problems were considered by |L2| Holtsmark (1919) 



in electrostatics and by || Chandrasekhar (1941) and PJTHI] Chandrasekhar & von Neumann 
(1942,1943) in a stellar context. We will show that many of the methods introduced by these 
authors can be extended to the case of point vortices, even if the calculations and the results 
differ due to the lower dimensionality of space (D = 2 instead of D = 3) and the different 
nature of the interactions. 

We consider a collection of N point vortices randomly distributed in a disk of radius R. 
We assume that the vortices have a Poisson distribution, i.e. their positions are independent 
and uniformly distributed over the entire domain. We are particularly interested in the 
"thermodynamical limit" in which the number of vortices and the size of the domain go to 
infinity (N — > oo, R — > oo) in such a way that the vortex density n = remains finite. In 



this limit, the Poisson distribution is shown to be stationary (see, e.g., [16] Novikov 1975) 
and is well-suited to the analysis of the fluctuations. For the moment, the vortices have the 
same circulation 7, but we shall indicate later how the results can be generalized to include 
a spectrum of circulations. 

The velocity V occurring at a given location of the flow is the sum of the velocities 
{i = 1, N) produced by the N vortices: 

N 

v = E*« (i) 

i=i 

*i = ~~r (2) 

where rj denotes the position of the i th vortex relative to the point under consideration 
and, by definition, rj_« is the vector rotated by +|-. Since the vortices are randomly 
distributed, the velocity V fluctuates. It is therefore of interest to study the statistics of 
these fluctuations, i.e. the probability W(V)d 2 ~V that V lies between V and V + dV . We 



1 The analogy between two-dimensional vortices and stellar systems is discussed in detail in 
Chavanis (1996,1998b) and § Chavanis et al. (1996). 
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find that this distribution behaves in a manner which is intermediate between Gaussian and 
Levy laws: the core of the distribution is Gaussian with "variance" 

(V 2 ) = ^\nN (3) 

while the high velocity tail decreases algebraically like V~ 4 . Since the "variance" behaves 
like ~ IniV, the thermodynamical limit is not well defined and the results are polluted by 
logarithmic corrections. Previous investigations of this problem were carried out numerically 
by @ Min et al. (1996) and @ Weiss et al. (1998) and analytically by Jimenez (1996). 



However, we must be aware that the knowledge of W(V) alone does not provide us with 
all the necessary information concerning the fluctuations of V. An important aspect of the 
problem concerns the speed of fluctuations, i.e. the typical duration T(V) of the velocity 
fluctuation V. This requires the knowledge of the bivariate probability W(V, A)d 2 Vd 2 A to 
measure simultaneously a velocity V with a rate of change 

dV N 

A =^T=E*< ( 4 ) 



* = -f(^-^^) (5) 



where v* = % is the velocity of vortex i. Then, the duration T{V) can be estimated by the 
formula: 

T(V) = -M= (6) 



where 



,„ / W(V, A)A 2 d 2 A 



is the mean square acceleration associated with a velocity fluctuation V. A similar quantity 
was introduced by P,|T0| Chandrasekhar & von Neumann (1942,1943) in a stellar context. 



We find that the distribution of the accelerations is governed by a Cauchy law and that the 
typical duration T(V) of a velocity fluctuation V behaves like V for V — > and V~ 1 for 
V — » oo. We also establish that the average duration of the velocity fluctuations is 

(T) ^= (8) 

ri7v In 

These results can be understood in the "nearest neighbor approximation" in which the most 
proximate vortex plays a prevalent role. In this approximation, the determination of the 
speed of fluctuations can be deduced from the theory of Smoluchowski (1916) concerning 
the mean lifetime of a state with X particles. 

In terms of the previous quantities, we can estimate the diffusion coefficient of point 
vortices by the formula: 
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D=-f T(y)W(V)V 2 d 2 V (9) 



We find that 



D~7v / hTiV (10) 



and we discuss qualitatively how the formation of "pairs" modifies the results of our study. 
In the context of freely decaying two-dimensional turbulence, the diffusion coefficient is time- 
dependent (since the circulation of a vortex increases as a result of successive mergings) and 
the diffusion is anomalous. From equation ([TUD, we establish a relationship 

"=i+| (ii) 

between the exponent of anomalous diffusion v and the exponent £ which characterizes the 
decay of the vortex density. This relation is in good agreement with laboratory experiments 
and numerical simulations. 

We indicate how our results are modified when we allow for a spectrum of circulations 
among the vortices. This is an important generalization since decaying 2D turbulence pos- 
sesses a continuous distribution of vortices. We show that the distribution of velocity and 
acceleration (or velocity gradients) are only slightly modified by the polydispersity of the 
vortices and we justify the validity of previous comparisons of full numerical simulations 
with vortex models that ignored this difference (e.g, [O] Jimenez 1996). 



Finally, we generalize our results to the case of vortex "blobs" with a finite core. We show 
that the natural cut-off at r = a, the vortex radius, removes the algebraic tail of the velocity 
distribution. Further, we show analytically that the distribution of accelerations makes a 
smooth transition from Cauchy (for small fluctuations) to Gaussian (for large fluctuations). 
It is likely that in between the distribution passes through an exponential tail as observed 
numerically by Min et al. (1996). 



II. THE STATISTICS OF VELOCITY FLUCTUATIONS 
A. The general expression for W(V) 

We shall now obtain the distribution Wn(V) of the velocity V produced by iV point 
vortices randomly distributed in a disk of radius R with uniform probability. To avoid a 
possible solid rotation, we shall assume that the system is "neutral" , the circulation of the 
vortices taking only two values +7 and —7 in equal proportion (the case of a non neutral 
system is treated in Appendix 0). Since a vortex with circulation —7 located in r produces 
the same velocity as a vortex with circulation +7 located in — r, and since the vortices are 
randomly distributed over the entire domain with uniform probability, the group of vortices 
with negative circulation is statistically equivalent to the group of vortices with positive 
circulation. We can therefore proceed as if there were a single species of particles but no 
solid rotation. Since we shall ultimately let R — > 00, we can assume without loss of generality 
that V is calculated at the center of the domain. 

Under these circumstances, the distribution Wn(V) can be expressed as: 
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N , N 



W N (V) = [ nr( ri )dV(v (12) 

J i=l ^ i=l ' 

where r(ri)d 2 Ti governs the probability of occurrence of the i-th point vortex at position r^. 
In writing this expression, we have assumed that the vortices are identical and uncorrelated. 
Now, using a method originally due to Markov, we express the 5-function appearing in flT2|) 
in terms of its Fourier transform 

*(x) = 77^ / e-**d 2 p (13) 



(2tt) 2 „ 

With this transformation, Wn(V) becomes 

W N (V) = ^J A N {p)e-^d 2 p (14) 



with: 



A N (p)=( f e^T^rj (15) 



where we have written 



* = -f% (16) 

If we now suppose that the vortices are uniformly distributed on average, then 

1 

vr'R 2 



= (17) 



and equation ( jl5|) reduces to 



Since 



7Ti? 2 V|r|=0 

we can rewrite our expression for A^ip) in the form 

N 



d 2 r = 1 (19) 



\ 7lR Z J\r\=0 J 



(20) 



We now consider the limit when the number of vortices and the size of the domain go to 
infinity in such a way that the density remains finite: 

N 

N —y oo, R — > oo, n = — — finite 

TlR 2 



If the integral occurring in equation (|20|) increases less rapidly than N, then 
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A(p) 



-nC(p) 



(21) 



with 



Cip) 



R 



|r|=0 



)d 2 r 



(22) 



We have dropped the subscript N to indicate that the limit iV —>■ oo, in the previous sense, 
has been taken. Note that A(p) can still depend on N through logarithmic factors, so that 
equation (|2T| ) must be considered as an equivalent of ( j2"0"| ) for large iV's, not a true limit. 

To calculate Cip) explicitly, it is more convenient to introduce $ as a variable of inte- 
gration instead of r. The Jacobian of the transformation {r} — > {$} is 



d(r) 



7 



47T 2 $ 4 



so that 



Cip) 



ry* r + OO 



or, alternatively, 



Cip) 



r 



4ir 2 J\®\ 



+oo 



;i-cos( P $))— d 2 $> 



Chosing polar coordinates with the x-axis in the direction of p, equation 
formed to 

r+oo r 2w 



ry" r+oo Hm rzn 



where 9 denotes the angle between p and Using the identity 



cos [z COS 



we obtain 



Cip) 



_ (1 - J (p*)) 



r 

2tt 7^ 

2ttH 
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or, writing x = p$, 



C(P) 



7 2 p 2 r+°° 



27T 



(1 - J (x)) 



dx 



(23) 



(24) 



(25) 

can be trans- 
(26) 

(27) 
(28) 



(29) 



Recall that Cip) must be evaluated in the limit N, R — > 00 with n = ^2 finite. Using the 
well-known expansion of the Bessel function Jq for small arguments: 



Jnix) 



x- 



1 - — + o ix* 
4 



(30) 
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we have the estimate: 



Cip) = ffU^L) (31) 
1D7T xwy^p 2 ' 



Since C(p) diverges weakly with N (logarithmically), the limiting process leading to formula 
(pl|) is permissible. For p > and N — > oo, we have 

A(p) = e-^ ln ^ 2 (32) 

and for p — ► 0, we obtain 

A{p) = e ^ lnp)p2 (33) 

The velocity distribution W(V) is simply the Fourier transform of A(p). We shall now derive 
the expression for W(V) in the core and in the tail of the distribution. 



B. The core of the distribution W(V) 

For V < V cr it(N), where V cr it(N) is defined by formula fl50D , the contribution of small 
p's in the integral (|i~4"|) is negligible and we can use the expression (|32|) for A(p). In that 
case, the distribution W(V) is the Gaussian 

W(V)= 2 \ e -^^ v2 (V<V cnt (N)) (34) 
ny z in JM 

If we were to extend this distribution for all values of V, we would conclude that its variance 

(V 2 ) = ^\nN (35) 

diverges logarithmically when iV — > oo. This result was noted by Jimenez (1996), |Tj| 
Min et al. (1996) and |2(J Weiss et al. (1998) who applied a generalized form of the Central 
Limit Theorem. In fact, the Central Limit Theorem is not strictly applicable here because 
the variance of the velocity created by a single vortex 

($ 2 )= / B r(r)$Vr= A* ^ J^jfcrrdr (36) 

diverges logarithmically when r — > 0; still, the distribution of V is Gaussian (for V % 
V cr u(N)) but its "variance" behaves like IniV. For V ^ V cr it(N), the distribution (j34|) 
breaks down because, for large velocities, the Fourier transform ( |i~4"D is dominated by the 
contribution of small p's and formula (|33| ) must be used instead of (|32|) . This implies that 
the high velocity tail of the distribution W(V) decays algebraically like V~ A (see section 
II C|) . This algebraic tail arises because we are on the frontier between Gaussian and Levy 
laws (see figure 1.1 of M Bouchaud & Georges 1990 and section |V|). 

The distribution (|34"D has been derived for a neutral system consisting in an equal number 
of vortices with circulation +7 and —7. If the system is non-neutral, there is a solid rotation 
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and the average velocity increases linearly with the distance. Therefore, at point a, equation 
([34]) must be replaced by (see Appendix |C]): 

W(V)= * e -^ {V -^ )2 (\V-± nia± \ <V crit (N)) (37) 
rrp In W 2 

The velocity distribution at a 7^ differs only from the distribution at the center of the 
domain by replacing the velocity V by the fluctuating velocity V = V — (V) = V — |n7a_|_. 
A factor 1/2 arises in front of the average vorticity nj because, for a solid rotation, the 
vorticity is twice the angular velocity. In Appendix ||, we extend our results to allow for a 
spectrum of circulations among the vortices. 



C. The high velocity tail of the distribution W(V) 

We shall now determine the behavior of the distribution W(V) for V — > 00. Introducing 
polar coordinates with the x-axis in the direction of V, and using (21), equation (|14D can 
be transformed to 

1 /*7r /"+oo 

W(V) = — / dO e- ipVcos9 e~ nC{p) pdp (38) 
2n 2 Jo Jo 

With the change of variables z = pV and t = — cos 9, equation (RSI) can be rewritten 

W(V) = -i— ^ e ! +1 -=£= [ + °° e^e- nC ^) zdz (39) 



2tt 2 V 2 J -1 ^r^PJo 

In this expression, t and z are real and the domains of integration tq : — 1 < t < 1 and 
Co : < z < +00 are taken along the real axis. Under these circumstances, the integral is 
not convergent if we expand the quantity 

in a power series of y, for V — > +00, and evaluate the integral term by term. However, 
regarding z and t as complex variables, it is possible to choose paths of integration along 
which this expansion will converge 0. 

We shall first carry out the integration on z, for a fixed t. It will therefore be possible 
to choose the (complex) integration paths for z dependent on t. The integration paths are 
modified as follows: To is replaced by r, the semi-circle with radius unity lying in the domain 
Im{t) > 0. Therefore, arg(t) varies from tc to when t moves from —1 to +1. On the other 
hand, Co is replaced by £ Wt , the line starting from the origin and forming an angle 

^ = Ki _arg(t) ) (40) 

with the real axis. When t moves from —1 to +1, u t varies from — ^ to For \z\ — > 00, 
according to (B3), we have: 



2 This method is inspired from || Chandrasekhar & von Neumann (1942). 
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-nC(f) = p ~^\nN- 



(41) 



Since the argument of z 2 is between — | and |, its real part is always positive and the 
convergence of equation d39|) is undisturbed. On the other hand, the argument of izt is 
equal to 



7T 1 

2 + - 



7T 



8 V2 



arg(t) ) + arg(t) 



and lies between || and ^f. Therefore, the real part of izt is always negative and the 
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function e decays exponentially to zero as \z\ — > oo. Therefore, with the new paths of 
integration r and ( wt , it is possible to expand the integrand of equation ( |39|) in power series 
of Tr, for V — > oo, and integrate term by term. When ^ 



v 



and we can write: 
W(V) = 



1 



2n 2 V 2 



7Z P 



dt 



'j Sir ' v^- 



0, we have, according to (|55|) ; 

(42) 



Azt 



1 + 



~8t7 



WF 2 



(43) 



Since this integral is convergent along any line on which the real part of izt is negative, we 
can replace the integration path by the line forming an angle 



71 



2 - ar §(*) 



(44) 



with the real axis. On this new integration path 

izt = —y, y real > 



and we get 



W(V) 



2n 2 V 2 



+i 



dt 



+oo 



57T 



In 



i vT^Jo 

iy\l y 2 n^ 2 Yn.t y 2 



VJt 2 V 2 8vr t 2 V 2 



+ ... 



e~ y x 

y 



t 2 



dy 



(45) 



where we recall that t is a complex variable and the integration has to be performed over 
the semi-circle of radius unity lying on the domain I m (t) > 0. Writing t = e ie , we find that 



+i 



dt 



i t 2 ^l^t 2 



0; 



+i 



dt 



i t*VT^¥ 



+i 



hit 



2tt 

dt = 

i t 4 v / r^ 3 



(46) 



Therefore, 



W(V) 



2Atx 2 V 4 Jo 

In this expression, we recognize the T-function 

r+oo 

T(n + 1) 



V dy 



- y y n dy 



(47) 



(48) 



with n = 3. Its value is T(4) = 6 and we finally obtain 

W{ y) = £^l (V>VUN)) (49) 

Therefore, for sufficiently large values of V, the velocity distribution W(V) decays alge- 
braically, like V~ 4 . In section IV), we give a physical interpretation of this result in terms of 
the nearest neighbor approximation. 

^From equations ([34]) and (|4"9"D , we can estimate the value of the velocity V cr u(N) for 
which the distribution W(V) departs from the Gaussian. This V cr u(N) is obtained by seeking 



the point where the two regimes (p4[) and (j49|) connect each other. Neglecting subdominant 
terms in IniV, one finds simply 

2 1 y*2 

V^N) ~ ln 1/2 (lniV) (50) 

This result shows that the convergence to a pure Gaussian distribution is extremely slow 
with N as emphasized by @ Jimenez (1996), flTJ Min et a/. (1996) and |(J Weiss et a/. 
(1998). Since the distribution W(V) decreases like \^~ 4 for V — >• oo, the variance of the 
velocity diverges logarithmically. 

Note, finally, that the distribution of V x , the x-component of the velocity is 

W(V X ) = -^==e~^^ (V < V crit (N)) (51) 

w (V') = rVE ( V ~ V crit( N )) (52) 
o7rl/° 



D. The formation of "pairs" 

The previous results should be all the more valid if the velocity V is calculated at a fixed 
point of the domain. In such a case, there is no restriction on the possible values of V since 
a vortex can approach this point with no limit producing extremely large velocities. The 
situation is different if V is now the velocity experienced by a "test" vortex. Indeed, if a 
"field" vortex approaches the test vortex below a certain distance, then a "pair" will form 
and our treatment, which ignores the correlations between vortices, will clearly break down. 
These pairs have been observed and studied numerically by [^UJ Weiss et al. (1998). 

We can simply estimate the typical separation below which a pair will form by comparing 
the velocity produced by a single vortex ^ with the typical velocity V typ = (^|MniV) 1//2 



produced by the field (see equation ^5|). This yields 



dcrit(N) = (nn In N)' 1 ' 2 (53) 



a distance slightly smaller than the interparticle distance by a factor ~ 1/yhxN . In the 
mathematical limit iV — > oo, there is no pair, since d cr n — > 0. This result is in agreement 
with the stationarity of the Poisson process when iV — > oo: if the vortices are initially 
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uncorrelated, they will remain uncorrelated. However, the convergence is extremely slow 
with N and close pairs will always form in realistic situations. As emphasized by |20 Weiss 



et al. (1998), a system of 10 3 — 10 5 vortices has a behavior which is a combination of both 
low-dimensional behavior, i.e. closed pairs, and high dimensional behavior described by 
traditional stochastic processes. 

A pair can be either a "binary" (rapidly rotating around its center of vorticity) when two 
vortices of the same sign are bound together, or a "dipole" (translating or rotating), when 
two vortices of opposite sign pair off. Of course, binaries and dipoles behave very differently. 
If the vortex is engaged in a binary with long correlation time, then, for practical purposes, 
the relevant velocity to consider is not its own velocity (which has a rotating component), 
but rather the velocity of the center of vorticity which is induced by the rest of the system. 
Therefore, a binary simply behaves like a single point vortex with larger circulation and 
relatively slow velocity (otherwise this means that the binary is itself engaged in a pair). It 
may be noted that, in the case of real vortices (with a finite core), the formation of binaries 
is replaced by merging events. By contrast, a dipole moves by itself and behaves like a kind 
of particle undergoing fast ballistic motion. Its velocity may be large but, since it creates 
a dipolar velocity field, the previous results cannot be applied directly and an appropriate 
treatment is required. 

We therefore expect that the velocity distribution (|34])(f49|) which ignores correlations 
between vortices will break down for V 3> V cr i t (N) since, in that case, the velocity is 
entirely due to the nearest neighbor and pairs form. In the following, we shall account for 
this failure by introducing a cuf-off at some V max , i.e. 

W(V) = (V> V max ) (54) 

This is the simplest modification that we can make to account for the formation of pairs at 
large velocities. For vortices with size a, the Gaussian distribution ( |54"[ ) is valid for all V's 
until the natural cut-off at V max ~ (see Appendix [D]). 



III. THE STATISTICS OF ACCELERATIONS 
A. The general formula for W(V, A) 

We are concerned here with the calculation of the bivariate probability Wn(V , A) to 
measure simultaneously a velocity V with a rate of change A = According to equations 
(|I|) and (^), V and A are the sum of N random variables <frj and ip i depending on the 
positions rj and velocities Vj of the point vortices. However, unlike material particles, 
the variables {rj, Vj}, for different z's, are not independent because the velocities of the 
vortices are determined by the configuration {rj} of the system as a whole. However, for our 
purpose, it is probably a reasonable approximation to neglect these correlations and treat 
{rj, Vj} (i = 1, N) as independent variables. We shall only describe qualitatively how the 
formation of pairs affects our results. 

When this decorrelation hypothesis is implemented, a straightforward generalization of 
the method used in section [II A| yields 



W 



N 



(V, A) = JL / e-^^A N (p,a)d^a (55) 
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with 



A N (p,a) = ( f R r°° r(r,v)e i( P*^d 2 rd 2 v) N (56) 

V /|r|=0 7|v|=0 / 



where we have defined 



$ = - — ^ (57) 



7 /v_l 2(r ■ v)r ± 



27r V r 2 r 4 



(58) 



and where r(r, v) denotes the probability that a vortex be in r with velocity v. According 
to our initial assumptions, the vortices are distributed uniformly on average and r(r) is 
given by equation (|T7|). On the other hand, their velocity distribution r(v) is given by 
equations ( j34D fl49|) of section [II]. However, due to the formation of pairs, this distribution 
must be modified at large velocities (see section [11 D| ). Instead of introducing a sharp cut-off 
t(v) = at v > v max , we shall assume for convenience that the Gaussian distribution fl34|) 
is valid for all velocities 0. Therefore, the probability that a vortex be in r with velocity v is 

1 4 _ 4 t „2 

r(r,v) = —- x - e (59) 

vrit^ In JM 

It is remarkable that the distribution ( [59] ) is formally equivalent to the Maxwell-Bolzmann 
statistics of material particles at equilibrium. Owing to this analogy, we can interpret the 
variance 

2 

^=^L i n N (60) 

47T 

as a kind of kinetic "temperature" . More generally, the moment of order p of the velocity is 
where the T-function is defined by equation ( fl8|) . In particular 

,2 



lnA^) 2 (62) 



V 16 

Recall that the distribution ( (5"9"D is valid only for a neutral system made of an equal 
number of vortices with circulation +7 and —7 (if the system is non neutral we must account 
for a solid rotation). Since a vortex with circulation —7 located in r and moving with 
velocity v produces the same velocity V and acceleration A as a vortex with circulation +7 



3 In Appendix [B|, we generalize our results for an arbitrary isotropic distribution of the velocity 
r(M). 
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located in — r and moving with velocity — v, and since the vortices are randomly distributed 
with uniform probability and isotropic velocity distribution, the two groups of vortices are 
statistically equivalent. Therefore, as in section [TT], we can proceed as if we had a single type 
of vortex with circulation 7 and no solid rotation. In Appendix ||, we extend our results to 
allow for a spectrum of circulations among the vortices (still for a neutral system) and in 
Appendix we consider the case of a non-neutral system. 

Substituting the expression for r(r, v) in equation fl55D, we obtain 



A N (p,a) = (-L l R f " -^r^e-^ v2 e^^d 2 vd\) N (63) 
\ttR 2 J\t\=o J\v\=o rry 2 In N J 

As in section 11 A , it is more convenient to use $ and "0 as variables of integration rather 
than r and v. The Jacobian of the transformation {r, v} — > {<I>.c'} is 



d(r, v) 7 4 



9($,#) 



16vr 4 $ 5 



(64) 



We must next express v — |v| in terms of our new variables $ and t*. According to equations 
(|57|) and (HD, we get 

* = a " d ^ = £3 («) 



Hence 



V 2vr$ 2 

Thus, in these new variables, the expression for A^{p,a) becomes 



(66) 



A »M = (^2 C IT -^e'^^^^TT^^^Y ( 67 ) 
\ttR 2 J^\=^ n Jpp\=o ny 2 In N 167r 4 $ 8 / 

The integral with respect to ip is Gaussian and is readily evaluated. We are left with 

A N (p,a) = (-L / + °° e^e-?- 1 "^ 2 * 4 -^^*^ (68) 

1 1 Zty ri 



We verify that 



1 /-+00 /y 



2 



<r* = 1 (69) 



7rR 2 y*| = _2_ 47T 2 $ 4 

Hence, we can rewrite equation ( |68| ) in the form 

1 1 ZTYJrC 

We shall now consider the limit N, R — > 00 with n = ^2 finite. If the integral appearing 
in equation (|70|) increases less rapidly than N, then 
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A(p,a) = e- nC( P& (71) 

with 

C(p,a) = ^- r°° (1 - e ^ e -|™ ln ^ 4 )i-d 2 $ (72) 

1 zirH. 



As in section [II A| , the function A(/?,o") represents an equivalent of Apf(p,o) for large N, not 
a true limit. It can therefore still depend on N through logarithmic factors. 

After introducing polar coordinates and integrating over the angular variable using (^T|) . 
we arrive at 

C(p,a) = f f °°(1 - J (p$)e-^ Mi )§ (73) 

Equations (|55|)(|7lD and (|73|) formally solve the problem but the integrals look difficult to 
calculate explicitly. However, if we are only interested in the moments of A for a given V 
(resp. in the moments of V for a given A), we only need the asymptotic expansion of A(p,a) 
for \a\ — > (resp. \p\ — > 0). This problem will be considered in sections [ill C| . First, we shall 
derive the unconditional distribution of the acceleration. 



B. The Cauchy distribution for A 

According to equation (p>5]), we clearly have 

W(A) = J-. / Afarfe-WWtfpdWV (74) 

l07T 4 J 

Using (]I3|), the foregoing expression for W(A) reduces to 

W{A) = ^J A{o)e-** A d 2 o (75) 

where we have written A{a) for A(0,<r). Hence, according to equations ([71]) and (0), we 
get 

A(a) = e- nC{a) (76) 

with 

^4r(l-e>«' -)* (77) 
Z7T Jo y? d 

Following the usual prescription, we have let R — > oo since the integral is convergent when 
$ — > 0. Integrating by parts, we obtain 

C(a) = ^-VnlnNa (78) 
8 

Hence 
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A(a) = e -h 2 " 3/2 ^° (79) 

The distribution W(A) is simply the Fourier transform of the exponential function ([T9"|). 
This is the 2D-Cauchy distribution: 

32 1 

W ( A ) = , n ,r (80) 

I 1 T n 3 7 4 lnA r I 

The Cauchy distribution is a particular Levy law. As such, it decays algebraically for large 
|A|'s, according to 



This result has a clear physical interpretation in the nearest neighbor approximation (see 
section [V]). Only the moments of order p < 1 of the acceleration exist and we have the 
general expression 

W ^(^^)^r(i^)r( 1 + |) m 



Note, finally, that the distribution of A x , the x-component of the acceleration, is the ordinary 
lD-Cauchy law 

W ^ = 3/2 2 A AT 1 i ~64A 2 ( 83 ) 

Clearly, the distribution of acceleration is related to the distribution of velocity gradients 
5V (see [[lj Jimenez 1996 and |5| Min ei a/. 1996) 



C. The moments (A 2 ) v and (^ 2 )a 

Let us introduce a Cartesian system of coordinates and denote by {A^} the different 
components of A relative to that frame. The average value of A^Ay for a given velocity V 
is defined by 

(A^A V ) v = / W(V, A)A,A v d 2 A (84) 

According to equation (|5"5]), it can be written 

WCVXVMv = ^-7 / A(p,a)e-^ V ^ A )A M A^V 2 ^ 2 A (85) 
or, equivalently, 

W(V)(^> V = ~J 4^{^ A) }WA (86) 
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Integrating by parts, we obtain 

1 r d 2 A 



W(V)(A„A U ) V = ~J ^^)e-^WA (87) 



Using the identity (|13|), we can readily carry out the integration on A and a to finally obtain 

W(V)(AMv = -±J ^^(p,0)e-^d 2 p (88) 

Since the characteristic function A(p,a) is isotropic in each of its variables [see equations 
([HD (|73|) 1, the tensor (A^A u )v is diagonal and can be expressed as 

(A,A u ) w = l -(A 2 ) w 5 tMV (89) 

where (A 2 )v is given by 

W(V)(A 2 )v = ~J T^y(A0)e^d 2 p (90) 

We need therefore the behavior of A(p,o) for \a\ — > 0, or according to (|7T|) the behavior of 
C{p,a) for \o\ — > 0. Introducing the function C(p) defined in section [11 A| [see equation (^8|)], 
and using equation (0), we can write 

C(p,a) = C(p)+F(p,a) (91) 

with 

F(p,a) = £- f°° J ( P m - e-frnW*)** (92) 

Z7T JO 3> 

We have let R — > oo since the integral is convergent when $ — > 0. In terms of this new 
function, the expression for (A 2 )v can be rewritten 

W(V)(A 2 ) V = £ / e- nC ^^(p, 0)e-^d 2 p (93) 

It is shown in Appendix [A| that 

f) fp 77-77 
^(p,0) = T7 2 lniV5(p) (94) 

where 5 stands for Dirac ^-function. Substituting the foregoing expression in equation (j93"D, 
we obtain 



47T 

|2 



iy(V)(A 2 ) v = ^lnAT (95) 



Therefore, (A }v behaves like the inverse of the velocity distribution W(V). Combining 



equations fl95|) (gj) and (f49|) , we obtain 
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(A% = ^(lniV)V^ (V <V mt (N)) (96) 

(A 2 ) v = nn\nNV 4 (V > V cnt (N)) (97) 

For V — > oo, (>l 2 )v behaves like l/ 4 . This result finds a simple physical interpretation in 
the nearest neighbor approximation (see section |V]). 

By a similar procedure, we can calculate the variance of the velocity for an assigned rate 
of change. We find: 

(^ 2 >A = ^lniV (A<A crit (N)) (98) 

OA 

(V 2 )a= . ] M (A>A crit (N)) (99) 
ixyn in JM 

For moderate values of A, the variance {V 2 )a is independent from A and coincides with 
formula (|35l). For large fluctuations, the result (|99| ) can be recovered in the nearest neighbor 
approximation (see section |V|). The crossover between the two distributions (|98| ) (|99| ) occurs 
at 

A crit (iV)~— ^(lniV) 3 / 2 (100) 
o 

a formula which also has a clear physical interpretation within the nearest neighbor approx- 
imation. 



IV. THE SPEED OF FLUCTUATIONS AND THE DIFFUSION COEFFICIENT 
A. The mean lifetime of a velocity fluctuation V 



On the basis of very general considerations (see, e.g., |2(J Weiss et al. 1998), we would 



expect that the typical duration of the velocity fluctuations be 



T typ ~ -jL= (101) 



This corresponds to the time needed by a vortex with typical velocity J (V 2 ) to cross the 
interparticle distance d ~ rr 1 ! 2 . Using expression (^) for (V 2 ), we find 

T typ )== (102) 

Apart from the logarithmic correction, this formula can be obtained by direct dimensional 
analysis. However, if we define the mean lifetime of a state V by the formula 
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the theory developed along the previous lines enables us to obtain a more precise characteri- 
zation of the speed of fluctuations depending on their intensity. This definition is consistent 
with the definition introduced by [0,rO] Chandrasekhar & von Neumann (1942,1943) in the 
context of stellar dynamics. Of course, equation (|103|) is just an order of magnitude but 
it should reasonably well account for the dependence on the speed of fluctuations with V. 
Using (|96|) and (|9~7D we find explicitly 

T (V) = J^l^ ^^ (V<V crit {N)) (104) 



T(V)= r =L= = 1 (V>V crit (N)) (105) 

For weak and large fluctuations, T(V) decreases to zero like V and V~ l respectively. These 
asymptotic behaviors are consistent with the theory developed by [18j Smoluchowski (1916) 
in his general investigation on fluctuation phenomena (see section |VC|) . These results (and 



those of section |TTJ) should be all the more valid if V is calculated at a fixed point. By 
contrast, if V denotes the velocity experienced by a test vortex, we expect some discrepancies 
at large V's due to the formation of pairs. In such a case, the correlation time can be 
extremely long (in particular for binaries). 

The average duration of the fluctuations is defined by 

r+oo 

(T) = / T(V)W(V)2itVdV (106) 
Jo 

To leading order in IniV (i.e. extending equations (|34]) and ( |104| ) to all V^'s), we obtain 

4/7T\ 1 / 2 1 

^ = 3V6J n 7v / hTiV (10?) 
in good agreement with the estimate ( |102j ) based on general physical grounds. 



B. The diffusion coefficient 

According to the previous discussion, we can characterize the fluctuations of the velocity 
of a point vortex (or a passive particle) by two functions: a function jy(V) which governs 
the occurrence of the velocity V and a function T(V) which determines the typical time 
during which the vortex moves with this velocity. Since the velocity fluctuates on a typical 
time T typ = ~/^f^ which is much smaller than the dynamical time Tjj = needed by 

the vortex to cross the entire domain, the motion of the vortex will be essentially stochastic. 
If we denote by P(r, t) the probability density that the particle be found in r at time t, then 
P(r, t) will satisfy the diffusion equation 
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dP 

— = DAP (108) 
dt v ; 



If the particle is at r = r at time t — 0, the solution of equation ( |108| ) is clearly 

J^(r,t[r ) = i _L_ e - l£ S^ (109) 

where D is the diffusion coefficient. The mean square displacement that the particle is 
expected to suffer during an interval of time At large with respect to the fluctuation time 
Ttypi is 

((Ar) 2 )=4DAt (110) 

We can obtain another expression for ((Ar) 2 ) in terms of the functions iy(V) and T(V) 
defined in the previous sections. Indeed, dividing the interval 

Ar= / V(t')dt' (111) 



into a succession of discrete increments in position with amount T(Vi) V«, we readily establish 
that 

((Ar) 2 ) = (T(V)V 2 )At (112) 



Combining equations (|110|) and (|112|) we obtain an alternative expression for the diffusion 
coefficient in the form 

D = 1 J T{V)W{V)V 2 d 2 V (113) 

Substituting for T( Vj and W(V) in the foregoing expression, we obtain to leading order in 
IniV: 

1 /6\V2 



D = -(-) 7 VmiV (114) 



We should not give too much credit to the numerical factor appearing in equation (|114|) 
since the definition ( |103|) of T(V) is just an order of magnitude. Note that the functional 
form of D is consistent with the expression 

D~T typ (V 2 ) ~7v / hTiV (115) 

that one would expect on general physical grounds. |2(| Weiss et al. (1998) have proposed 
to describe the motion of the vortices by a Ornstein-Uhlenbeck (OU) process with "friction" 
— V/ < T > to take into account the finite decorrelation time of the system. The function 
T(V) introduced in the present article could be used instead of (T) to take into account the 
dependence of the decorrelation time with the strength V of the fluctuations. 

The present theory ignores the formation of pairs since we have formally introduced a 
cut-off at large Vs. This cut-off is justified for binaries since, as we have already explained, 
the relevant velocity to consider is the velocity of the center of vorticity not the velocity of 
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the individual vortices engaged in the pair. This is not the case for dipoles which can make 
relatively long jumps from one point to another with an almost ballistic motion. |HJ Weiss 
et al. (1998) have proposed to interpret these jumps in terms of Levy walks responsible 
for anomalous diffusion. Therefore, the present theory should provide reliable results only 
for relatively short times. For times t ^> To, the pairs must be taken into account and 
anomalous diffusion will ensue. The case of passive particles is not so different. A passive 
particle is advected by the other vortices but has no influence on their motion. However, 
passives can get trapped in the vicinity of a vortex undergoing fast dipolar motion and 
experience Levy walks also. 

In |T7| , we have actually made the connection between Levy walks and anomalous diffu- 
sion more explicit, by showing that the flight time p.d.f. has a slow power-law decaying tail 
with an exponent that can be related to the number of vortices decay exponent £ (see next 
section). 



C. Application to 2D decaying turbulence 

Let us consider a collection of vortices of size a, vorticity to and density n. Due to merging 
events, their size increases with time as the density decreases. The typical core vorticity 
uj remains constant during the course of the evolution as suggested in H Carnevale et al. 



(1991) [|TJ Benzi et al. (1992) and |L9| Weiss et al. (1993), and observed experimentally 
by [T1J Hansen et al. (1998). These authors found that the density decreases as n ~ t~^ 
with £ m 0.7. As the energy E ~ Nu 2 a A is conserved throughout the merging process, the 
typical vortex radius a ~ Since the average distance between vortices, of order d ~ t^ 2 
increases more rapidly than their size, the point vortex model should provide increasing 
accuracy. We can therefore expect that the vortices will diffuse with a coefficient D ~ 7 



see equation (|114 )1 where 7 ~ ua is their circulation (we ignore logarithmic corrections). 



If the diffusion coefficient were constant, then the dispersion of the vortices 

(r 2 ) ~ Dt (116) 

would increase linearly with time as in ordinary Brownian motion. However, since D varies 
with time according to 

D ~ Lua 2 ~ t^ 2 (117) 

we expect anomalous diffusion, i.e. 

(r 2 )-^ (118) 



with v ^ 1. Substituting equation (|117 ) into equation ( |116|) we obtain the following relation 
between v and £: 

v=\ + \ (119) 



This expression differs from formula (19) of JTTJ] Hansen et al. (1998) because their 
estimate of D is different. These authors estimate the diffusion coefficient by D ~ T merg (V 2 ) 
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where T merg is the average time between two successive mergings of a given vortex. By a 
simple cross section argument, they obtain r mero ~ — y= which is larger than Tn uct ~ 

a ny/ (V 2 )a J 

— by a factor -. This shows that the merging time is not the relevant correlation time 

■J(v 2 ) a 



to consider in the diffusion process. In fact, using £ m 0.7, formula ( |119| ) leads to v « 1.35 in 
better agreement with the experimental value v m 1.3 ( v w 1.4 for passive particles) than 
their relation v — 1 + ^ which leads to v « 1.53. Formula ( |119| ) is also in perfect agreement 
with the numerical simulations of |T7] Sire & Chavanis (1999). 



V. THE NEAREST NEIGHBOR APPROXIMATION 

A. The importance of the nearest neighbor 

The velocity V and the acceleration A experienced by a test vortex (or occurring at 
a fixed point) are the sum of iV random variables and ^ produced by all the vortices 
present in the system [see equations (0) and (|)]. In each sum, the highest term is due to 
the nearest neighbor, at an average distance d ~ n -1 / 2 from the point under consideration. 
This single vortex creates a typical velocity and acceleration 

VL ~ fiV ~ £z£ (120) 



.2nd J Att 2 7tR 2 



A 2 ~SL_^(JL) (121) 

\2tt<P) 4vr 2 URV { ' 

It is interesting to compare the contribution of the nearest neighbor with the contribution 
of the other N — 1 vortices. For that purpose, we estimate the typical value of V and A 
produced by all the vortices by 



r 2 . . ATI * \ _ . AT I ' „?2_ I 



\4 7r 2 r 2/ _/ |r|=d y 47T 2 r 2 47T7Ti? 2 

A 2 ~ Ntf(-X-r) ~ ^ r(r)-^-rf 2 r ~ £^(JLV (123) 
\47r 2 r 4 / i[ r |=d v ; 4?r 2 r 4 4tt Vvri? 2 / v ; 

If the variance of $ and "0 were finite, the Central Limit Theorem would be applicable 
and the variables V and A would scale like \/~N. In that case, none of the terms in the 
sums ([[]) (§J) would have a dominant contribution and the scaling a/ZV would simply reflect 
the collective behavior of the system. This is not the case however in the present situation 
since the variance of <& and tp diverge. The variance of if) diverges algebraically and, thus, 
the acceleration produced by all the vortices is dominated by the contribution of the nearest 
neighbor. This individual nature is a specific and striking property of a Levy law. For a Levy 
law, the sum of N random variables behaves like the largest term (compare equations ( |1 23| ) 
and (|121| )). The case of the velocity is particular because the variance of diverges only 
logarithmically. First considering equation ( |122| ), and neglecting the logarithmic correction, 
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we observe that the velocity produced by all the vortices behaves like v~N, as though the 
Central Limit Theorem were applicable. However, comparing with equation ( |120| ), we note 
that \/N is also the scaling of the largest term in the sum. Therefore, the velocity has a 
behavior which is intermediate between Gaussian and Levy laws, as remarked earlier. In 
a sense, we can consider that the velocity is dominated by the contribution of the nearest 
neighbor and that collective effects are responsible for logarithmic corrections. 



B. The distribution due to the nearest neighbor 

In light of the previous discussion, it is interesting to analyze in more detail the distri- 
bution of the velocity and acceleration produced by the nearest neighbor. For that purpose, 
we must first determine the probability T nn (r)dr that the position of the nearest neighbor 
occurs between r and r + dr. Clearly, r n . n (r)dr is equal to the probability that no vortices 
exist interior to r times the probability that a vortex (any) exists in the annulus between r 
and r + dr. Therefore, it must satisfy an equation of the form: 

r n .n(r)dr = (l - £ T n . n {r')dr^jn2nrdr (124) 



where n = denote the mean density of vortices in the disk. Differentiating with respect 
to r we obtain 



d 
dr 



-r n . n (r) (125) 



2ixnr 

This equation is readily integrated with the condition r n n (r) ~ 2imr as r — > 0, and we find 

T n . n (r) = 2Trnre- nnr2 (126) 

This is the distribution of the nearest neighbor in a random distribution of particles. From 
this formula, we can obtain the exact value for the "average distance" d between vortices. 
By definition, 

d= / r n . n (r)rdr (127) 
Jo 

Hence 

d=^-= (128) 
2v^ 

The probability of finding the nearest neighbor in r with velocity v is 

r„. n (r, v) = ne"™- 2 x i e ~^W v2 (129) 

n^y z In TV 

If we assume that the velocity V and acceleration A are entirely due to the nearest neighbor, 
then 

W n . n (V, A)d 2 Vd 2 A = r n . n (r, v)rf 2 rd 2 v (130) 



22 



with 



V = "££ (131) 
A=-^fe-^i) (132) 



27T V r 2 

Since the Jacobian of the transformation {r, v} — > {V, A} is 

d(r, v) 7 4 



a(v,A) 



16tt 4 V 8 



(133) 



we obtain: 



where we have used 



^ 2 1 «7 2 A 2 

WUY > A) = I^hTiV yse-^e-^** ( 134 ) 



and v = TTTF^ ( 135 ) 



2ttK 2tt \/ 2 

The nearest neighbor approximation is expected to give relevant results only for large values 
of the velocity and the acceleration. Thus, we can make the additional approximation 

^.»(V,A) = i ^i 5 e-^ (136) 
Integrating on the acceleration, we find 

W-CV) = T^yi (137) 



in perfect agreement with equation PS|) valid for 1/ > V^. This shows that the algebraic 
tail of the velocity distribution is produced by the nearest neighbor. This is characteristic 
of a Levy law. On the other hand, for V < V cr u, the velocity distribution is Gaussian as 
if the Central Limit Theorem were applicable. Once again, the simultaneous occurrence 
of collective and individual behaviors is a manifestation of the very peculiar nature of an 
interaction in r~ l in two dimensions. 
Integrating on the velocity, we find 

WUA) = (138) 



in perfect agreement with the asymptotic behavior of the Cauchy distribution (|81|) . We also 
establish that 

(A 2 ) v = rwrln NV 4 (139) 
2 A 

(V 2 )a = —== (140) 
7T\/nmi\ 



in complete agreement with formulae (971) and (199]). 
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C. The application of Smoluchowski (1916) theory 



In the nearest neighbor approximation, the duration of the velocity fluctuations can be 
deduced from the theory of fl8| Smoluchowski (1916) concerning the persistence of fluctu- 
ations. This approach was used by H Chandrasekhar (1941) in his elementary analysis of 
the fluctuations of the gravitational field. An account of Smoluchowski theory can be found 
II | Kandrup (1980). In the case of point vortices, it leads to the formula 



in 



7 V 

T{V) = g- (141) 

v 1 4t> vdL + v 2 

where V = is the velocity due to the most proximate vortex, at a distance r from the point 
under consideration. The Smoluchowski formula (|141|) has the same asymptotic behaviors as 
( |104j ) and ( |105| ). These asymptotic behaviors have a clear physical meaning. When r = ^7 
is small, corresponding to large velocities, it is highly improbable that another vortex will 
enter the disk of radius r before long. By contrast, on a short time scale T ~ = ~ ^7, the 
vortex will have left the disk. When r = 7^7 is large, corresponding to small velocities, 
the probability that the vortex will remain alone in the disk is low. The characteristic 
time before another vortex enters the disk varies like the inverse of the number of vortices 
expected to be present in the disk, i.e. T ~ =—^r ~ -^=. The demarcation between weak 

and strong fluctuations corresponds to V ~ ^n 1 / 2 , i.e. to the velocity produced by a vortex 
distant rr 1 ! 2 from the point under consideration. 



VI. CONCLUSION 

In this paper, we have analyzed in detail the statistics of velocity fluctuations produced 
by a random distribution of point vortices. We have determined the velocity distribution and 
the speed of fluctuations. We have also shown how some of the results can be understood 
in the nearest neighbor approximation. In a future study (in preparation), we shall be 
concerned with the spatial and temporal correlations of the velocity. Our results should be 
accurate if the velocity is calculated at a fixed point. However, there should be substantial 
discrepancies at large velocities if V now represents the velocity experienced by a point 
vortex. This is due to the formation of pairs (binaries or dipoles) when two vortices come 
into contact. If we ignore these pairs, the motion of the vortices is purely diffusive and we 
determined the functional form of the diffusion coefficient (up to a numerical factor). In 
the case of real vortices, with a finite core, the formation of binaries is replaced by merging 
events and the number of vortices decreases with time. This results in anomalous diffusion. 
We proposed a relationship between the exponent of anomalous diffusion v and the exponent 
£ which characterizes the decay of the vortex density. 

Our study was directly inspired by the work of H,|K| Chandrasekhar & von Neumann 
(1942,1943) concerning the statistics of the gravitational field arising from a random distri- 
bution of stars. Of course, the results differ due to the different nature of the interactions 
but the mathematical structure of the problem is the same. This is another example of the 
analogy between two-dimensional vortices and stellar systems investigated by [^,0 Chavanis 
(1996,1998b) and |7j Chavanis et al. (1996). Likewise in the stellar context, the formation 
of binaries alters the results of the stochastic analysis at large field strengths. 
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APPENDIX A: DERIVATION OF FORMULA (91). 

By definition, 

F{p,a) = f r J (p*)(l - e- l ^ N ^)§ (Al) 

Z7T JO <P 

For p = 0, we have already found [see equations ( [77|) and ([7^)1 that 

7 2 / 

F(0,a) = —Vn In No (A2) 
8 

Therefore: 

(p, 0) = +oo if p = (A3) 



d(a 2 , 

For p ^ 0, we can make the change of variable z = p$ in equation ( |A1|) . This yields 

F^) = ^/ + ^ oW(1 _ e -J— -J*,* (A4) 

We have therefore to determine the behavior of the function 

f(p)= [ + °° Jo(z)(l - (A5) 
jo z a 

as p — > 0. Clearly, it is not possible to expand the quantity 

1 - e~ pzi 

which occurs under the integral sign as a power series of pz 4 and evaluate the integral term 
by term. However, writing the Bessel function in the form 

J °w-M>7rb (A6) 

and regarding z and t as complex variables, it is possible to choose integration paths along 
which this expansion will converge. Using the contours introduced in section |II Q , the 
function f(p) can be rewritten 

We readily verify that the real parts of izt and —pz 4 are always negative, so the convergence 
of equation ( |A"5| ) is not disturbed. With these new contours, it is now possible to expand 
the integrand in a power series of pz 4 and integrate term by term. For our purposes, it is 
only necessary to consider the term of first order in this expansion 

f(p) = p-K e f -=£= / e M zdz + 0(p 2 ) (A8) 
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Therefore 

/'(0) = -K e f -y£= I e^zdz (A9) 

The integration on z can be carried out equivalently along the line (^ t , defined in section 
II C , on which izt = —y, y > real. We obtain 

f(0) = --n e -7== 2 e-y^dy (A10) 

where t is a complex variable running along the semi-circle of radius unity lying in the 
domain T m (t) > 0. Since 

; = All 



we find /'(0) = 0. Therefore 

<9F 



d(a 2 , 

To prove formula (05), it remains to show that 



(p,0) = if p^O (A12) 



0F >,0)d 2 p=^ 7 2 lniV (A13) 



For that purpose, we introduce the function 

/(a 2 ) = J F(p,a)d 2 p (A14) 
Substituting explicitly for F(p,a) and introducing polar coordinates, we get 

r+oo r+oo , _ . ^/<J) 

/(a 2 ) = 7 2 / prfp / J (p$)(l - e-l™ 1 *^ 2 * 4 )^ (A15) 

Under this form, it is not possible to interchange the order of integration. An alternative 
expression for /(cr 2 ) can be obtained along the following lines. Writing 

" + 0O 



where 



r+oo r+oo 

I(a 2 ) = 7 2 / dp p$J (p$)g($)d$ (A16) 
Jo jo 



g(9) = (1 - e -^ N ^ 4 )^ (A17) 



and integrating by parts with the identity 



xJo(x) = -^{xJi{x)) (A18) 
dx 



we obtain 
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f+OO f + OO 

Ua 2 ) = - 7 2 / dp / Ji(p$)$</'($)d$ (A19) 

JO JO 



It is now possible to interchange the order of integration. Since 

"+OO 



f+OO 

/ J x {x)dx = 1 (A20) 

JO 

equation ( |A19| ) reduces to 

r+oo 

/(a 2 ) = - 7 2 / g'(1>)d* = 7 2 g(0) (A21) 

JO 



Hence 



7(^2) = ^flniVa 2 (A22) 
This formula is valid for any value of a but, for our purposes, we only need the result 

7'(0) = ^ 7 2 lnAT (A23) 

Since, by definition, 

OF 



m = Jd^^ d2 P (A24) 
we have proved equation ( |A13|) . 

APPENDIX B: GENERALIZATION TO INCLUDE A SPECTRUM OF 
CIRCULATIONS AND AN ARBITRARY ISOTROPIC DISTRIBUTION OF 

VELOCITIES 

So far, we have assumed that the system was a "vortex plasma" consisting of an equal 
number of vortices with circulation +7 and —7. We shall now indicate how the previous 
results can be extended to include a spectrum of circulations among the vortices. We 
shall also relax the assumption (|59| ) concerning the velocity distribution of the vortices 
and generalize the results of section |TXT| to any isotropic distribution r(v) = t(|v|) of the 
velocities. Such a distribution can be written conveniently in the form 



t(v) = / ^S(v - v )d\ (Bl) 

J ZTCVq 

If r(7) governs the distribution over the circulations, it is clear that equation (K^) has 
to be modified according to 

A N (p,o) = ( [ +CO [ R / + °° r(7)r(r)r(v)e l ^ + °^)rf7rf 2 rc/ 2 v") JV (B2) 

\J7=-oo J|r|=0 J|v|=0 / 

There is no a priori restriction on the function r(7) but we shall be particularly interested 
in the case where the system is "neutral", i.e. 7 = / r^jdj = 0. It is only in this 
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circumstance that the velocity distribution ( [BID ma y De used. Otherwise, there is a solid 
rotation of the system which adds to the dispersion of the particles (see Appendix |C|). 
The expression ( pTj) for A N (p,a) is now replaced by 

/ ] r+oo r+oo 

A N (p,a) 



K nR 2 y 7 =-^ y|v i=o ii*i=^L y^| =0 

*^ T MiB - y^^^^y <B3) 



Introducing polar coordinates, using identity (|27|) and substituting for 

\j\if) \ 2tt$ 2 / 2tt$ 2 ' 



5 



(sW-^-TTT-^Tr) (B4) 



17 

in equation (|B3|), we can easily integrate on ip to obtain 



/ 1 r+oo r+oo r+oo f^ixn \ '"V 2 



JV 



(B5) 



r+oo r+oo r+oo 

/ / L r (7) r ( v o)^^7 d~fd\ d 2 <f> = 1 (B6) 



It is readily verified that 

1 r+oo r+oo r+oo ^,2 
7lR 2 Jj=-oo J\vo\=0 M\=£L 

Therefore, the expression for A^{p,a) can be rewritten equivalently 

/ ] r+oo r+oo r+oo 

A N (p,a) = 1 - — / / / 

V TcR 2 J 7 =-oo 7|v|=0 J|< 



|v|=o7|$|=^L 

xr( 7 )r(v)(l- e ^J (^$ 2 )) i -^ I d 2 ^ 2 vd 7 )' V (B7) 



where we have written v instead of vo as it is a dummy variable of integration. In the limit 
N, R — > oo with n = ^2 fixed, we obtain 



A(p,a) = e - nC( fi' a) (B8) 

with 

For a = 0, the function C(p,a) reduces to: 



2 2 

C(p) In N (p>0) (BIO) 

107T 



C(p)~-^lnp (p-0) (Bll) 
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Therefore, the velocity distribution becomes 



W(V)= 4 e (V<V cnt (N)) (B12) 

nj 2 in JM 



^(V)~ (V>V mt (N)) (B13) 

with 

/n Z2 sl/2 

V crit (N) ~ f^lniVj (B14) 

These results differ from (|3~4D (f49[) (p0|) simply by the substitution 7 2 — > 7 2 where 7 2 = 
/ t(7)7 2 ^7 is the average enstrophy. 
Writing 

C(p,a)=C(p)+F(p,a) (B15) 

we get 

r+oo r+oo r+oo f fllXO \\ '"Y 2 

F(p,<r) = / r( 7 )r(v)J (p$) 1 - J ( -ppt/* 2 ^ ^ 2 vd$ (B16) 

J 7 =_oo i|v|=o ^0 V \ |7| //ZTr'P' 3 

Following the same steps as in Appendix [A], we find that 

f)p 

P 0) = -n\H{p) (B17) 



d(a 2 ) 
where 



w 2 = / T (\v\)v 2 2nvdv (B18) 

JO 

is the mean square velocity of the vortices for the isotropic distribution t(|v|). Therefore, 
equations (|96|) and (|97| ) are modified according to 

n 7 Z 



(^ 2 ) v = ^^^iniVe^in^ (y < ^(iV)) (B19) 



and equations 



T(V) 



T 

according to 
2V 



(B20) 



n y 7 



v/hTiVA 



(B21) 
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T(V) = - 

2nVv 2 V 

The mean duration of the fluctuations is 



(V>V crit (N)) 



(T) 



3v^6 V nv 2 



and the diffusion coefficient 



D 



1_/6\V 2 



72 V 



7T 



For a Gaussian distribution of the velocities, we recover the results of section |TTT 
ately generalized to account for a distribution over the circulations 
For p= 0, the function C(p,a) reduces to 

"+oo r+oo r+oo 



C(a) 



r( 7 )r(v)(l-J ( 



' 7=— oo J | v|=0 ■/ 

Integrating by parts and using the identity 

"+°° JAx 



2na ~ 2 
-r-rV® 

7 



■dx = 1 



we obtain 



where 



o x 



r 



(|v|)f 2irvdv 



is the average velocity of the vortices. Equation (|75 ) is changed to 

71. W I — 

A(<r) = g~ 2" 

and the Cauchy distribution 



to 



W(A) = =^ 



1 



7rn 2 |7| u 2 i 4A 2 

\ |7| f 



3/2 



We find also 



(V 2 )a 



(V 2 )a 



47T 

7* A 



27r| 7 |t> 



IniV (A < AritW) 
(A> A crii (iV)) 



with 



Acrit(N) ~ ~| 7 |mJln./V 



(B22) 

(B23) 

(B24) 
appropri- 

(B25) 

(B26) 
(B27) 
(B28) 

(B29) 
(B30) 



For a Gaussian distribution of the velocities, we recover the results of section pi 
ately generalized to account for a distribution over the circulations. 



(B31) 
(B32) 

(B33) 
appropri- 
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APPENDIX C: THE CASE OF A NON NEUTRAL SYSTEM 



In this appendix, we consider the case where all vortices have the same circulation 7 and 
we allow properly for the solid rotation of the system that this distribution involves. We 
shall be particularly concerned with the expression of the bivariate distribution W^V, A) at 
a point a rotating around the center of the domain with velocity (V)(a) = jn-ya±. 

It is clear that equations (|55|) (|56f) ([57]) and (|Blj|) remain valid in this more general 
situation if we make the substitutions r — > £ = r — a (relative distance) and v — > v' = 
v— (V)(a) (relative velocity). According to equation ( j37j ) that we shall prove ulteriorly, the 
probability that a point vortex be found in r with velocity v is: 

r(r, v) = -L x *^r^t^^f (CI) 

We find it convenient to introduce the variable 

w = v - (V)(r) = V + (V)(a) - (V)(r) = v' - l -n^ (C2) 
In terms of this new variable, the acceleration (^) becomes 

€-v'! 0(6 w) ■ - 2 >r^. (C3) 
We can restore the initial form of A(p,<?) by defining p 1 = p — \n^a^ in (|56|). Taking p' as 



a variable of integration in (|55|), we find that the original system (|55|) ( |56D fl5?|) (^8|) ( |59| ) 
is preserved if we make the substitutions r — > ^, v — > w and A — > A — 7^7 (in the 
following, the dummy variable is renamed p) . 

The characteristic function A(p,a) can be written 

A(p,a) = e - nC « < PW (C4) 

with 

Cafea-) = / r(w) (l - e'«)^w (C5) 

Since we are not at the center of the domain, the imaginary part of C a (p,a) does not vanish. 
We have: 



C a (p,a) = C(p,a)-zC'( Pl a) (C6) 



with 



C(fi,a) = J r(w)(l - cos(p$ + o#))d 2 Crf 2 w (C7) 

( "(p.a) = y r( w) sin(p5» + <nl»<f 2 Z<! 2 w (C8) 

The integrals C(p,a) and C'(p,a) have a very different structure. The first integral is domi- 
nated by relatively proximate vortices and can be approximated with sufficient accuracy by 
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the function evaluated in section pT[ In the second integral, the contribution of proximate 
vortices vanishes by symmetry. As a result, the integral is dominated by large values of 
|£| or, equivalently, small values of |#| and . We can therefore make the approximation 
sin(/ofl> + mp) ~~ (fb ai: and write 

C a (p,a) = C(p,a) -ipj $d 2 £ -ia J r(w)#d 2 &/ 2 w (C9) 

The last integral cancels out due to the isotropy of the velocity distribution r(w). In the 
second integral, we recognize the mean-field velocity created in a by the average distribution 
n of the vortices: 

(V)(a) = n f $d 2 £ = -n 7 a ± (CIO) 



Therefore: 

A(p,a) = e -"C(p.fr) + i P (v)( a) (cn) 
Combining our previous results, we obtain 

W(V, A) = -i- / e~ t( P V ^A(p,a)d 2 pd 2 <r (C12) 

167T 4 J 

where A(p,o) is the function considered in section JTT[ The bivariate distribution flC12p at 
differs only from the distribution at the center of the domain by replacing the velocity 
V by the fluctuating velocity V = V — (V)(a) = V — ^n-y&± and the acceleration A by the 
fluctuating acceleration A. = A — ^n r y\±. We can account for a spectrum of circulations 
among the vortices simply by changing 7 in 7 in the new terms and using the results of 
Appendix ||. 

The mean life of a velocity fluctuations V can be defined by 

T(V) = - 7 J2L= (CIS) 

and the diffusion coefficient by 

D = ^J T(V)W(V)V 2 d 2 V (C14) 



They have the the same functional form as (|104j)(|105|) and (|114 ). Therefore, the results are 



unchanged in the case of a solid rotation. For a differential rotation, the fluctuation time is 
related to the local shear S as investigated by || Chavanis (1998a) using an approximation 
in which the point vortices are simply transported by the mean flow. Therefore, the present 
theory gives the value of the fluctuation time in regions where the shear cancels out. Of 
course, a more general calculation should take into account the effect of both the shear and 
the dispersion of particles, but this will not be considered in that paper. Note also that 
when the system is inhomogeneous, the point vortices are subjected to a systematic drift ( 
Chavanis 1998a) in addition to their diffusive motion. This drift can be understood as 
the result of a polarization process. 
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APPENDIX D: THE CASE OF VORTEX BLOBS 



In reality, the vortices have a finite radius a which is not necessarily small. This finiteness 
is responsible for a maximum allowable velocity V max ~ achieved when two vortices are 
at distance ~ 2a from each other. Higher velocities are forbidden because the subsequent 
evolution is marked by merging events (see, e.g., Sire & Chavanis 1999). It is interesting to 
consider the distribution of velocity V and acceleration A produced by a collection of vortex 
"blobs" with finite radius. This problem was previously treated by Jimenez (1996) and 
|15|j Min et al. (1996) using numerical methods. The theory developed in this article makes 



possible to obtain new analytical results. 

Introducing a cut-off at r = a, equation (E3) is replaced by 



^w^/** 1 < D1 > 



2ttR 



In the limit N, R — > oo with n = -J*! finite, we obtain 

GW = ^! ln («) (D2 ) 
87r V a J 

In the case of vortex blobs, the singularity at p = is removed and the characteristic function 
Cn(p) is exactly quadratic. Therefore W(V) is the Gaussian ([34]) for all V < V max . There is 
no algebraic tail in the limit considered. However, the convergence to the limit distribution 



([34]) is still slow (see |l^] Jimenez 1996) and, in practice, the velocity distribution can differ 
substantially from the Gaussian even for large values of TV (note that the formalism presented 
in this paper could be used to study the dependence of the velocity distribution with the 
number N of vortices). 

More interesting is the distribution of the acceleration A. In the case of vortex blobs, 
equation (|77|) is replaced by 

C (a) = f- f^(l- e 4™^ 2 * 4 )|?j (D3) 



2ix Jo $ 3 
After integrating by parts, one obtains 



C(a) = vra 2 fe~lS^ ff2 - l) + -L Vn In iVErf Vnn In No) a (D4) 



Erf(x) = -^= f e~ y2 dy (D5) 

V 71 JO 



where 



is the "Error function" . 

For a — > 0, the function C(p) is quadratic in a: 

„. , n7 4 lnA r cr 2 , , .„ . 

d'^ ' (D6) 

implying that the distribution W(A) is Gaussian for large values of the acceleration: 
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W(A) ~ " A e"^^ A (|A| -> +00) (D7) 
n^7 4 IniV 



Its variance is: 



n 2 7 4 



167r 2 a 2 

For a — > +00, the function C(a) is linear in cr: 

7 2 



(A 2 ) = IniV (D8) 



C(a) ~ -WnlniVa (|a| -> +00) (D9) 



and we recover the Cauchy distribution ( PUD for small values of |A|. Therefore, the distri- 
bution W(A) makes a smooth transition from Cauchy (concave in a semi-log plot) for small 
fluctuations to Gaussian (convex in a semi-log plot) for large fluctuations. It is likely that 
in between the distribution passes through an exponential tail as observed numerically by 
[15| Min et al. (1996). Of course, when a is reduced, the transition between the two regimes 



happens at larger fluctuations (see equation ( |D8| )). According to equations fl52| ) and 
the relevant non-dimensional parameter to consider is the area fraction na 2 . In decaying 
turbulence, the influence of an extended core should be manifest at the beginning of the 
evolution. 
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